arXiv:1503.01899v2 [cond-mat.mes-hall] llJun2015 


Spin-Motive Forces and Current-Induced Torques in Ferromagnets 

Kjetil M. D. Hals 1 and Arne Brataas 2 
1 Niels Bohr International Academy and the Center for Quantum Devices, 

Niels Bohr Institute, University of Copenhagen, 2100 Copenhagen, Denmark. 

2 Department of Physics, Norwegian University of Science and Technology, NO -7 f 91, Trondheim, Norway. 

In metallic ferromagnets, the spin-transfer torque and spin-motive force are known to exhibit a 
reciprocal relationship. Recent experiments on ferromagnets with strong spin-orbit coupling have re¬ 
vealed a rich complexity in the interaction between itinerant charge carriers and magnetization, but 
a full understanding of this coupled dynamics is lacking. Here, we develop a general phenomenology 
of the two reciprocal processes of charge pumping by spin-motive forces and current-driven mag¬ 
netization dynamics. The formalism is valid for spin-orbit coupling of any strength and presents a 
systematic scheme for deriving all possible torque and charge-pumping terms that obey the sym¬ 
metry requirements imposed by the point group of the system. We demonstrate how the different 
charge pumping and torque contributions are connected via the Onsager reciprocal relations. The 
formalism is applied to two important classes of systems: isotropic ferromagnets with non-uniform 
magnetization and homogeneous ferromagnets described by the point group C 2 v 


I. INTRODUCTION 

A current of electrons that enters a ferromagnet can 
produce a torque on the magnetization due to the trans¬ 
fer of spin angular momentum from the itinerant car¬ 
riers to the magnetic system. Following the first theo¬ 
retical predictions reported in 1996 by Berger and Slon- 
czewskipB] the spin-transfer torque (STT) has garnered 
abundant experimental evidenced Today, current-driven 
magnetization dynamics is a highly active research area 
with a large application potential in magnetoelectronic 
devices (e.g., the first STT-RAM entered the market¬ 
place in 2012). EO 

The current-driven dynamics of metallic ferromagnets 
is described by the Laudau-Lifshitz-Gilbert-Slonczewski 
(LLGS) equation : 3 

m = — y m x H 0 ff + m x rmi + T. (1) 

Here, m(r, t) is a unit vector along the magnetiza¬ 
tion M(r,<) = M s m(r,f), and 7 is the gyromagnetic 
ratio. The first term on the right side of Eq. ([T]) 
is the torque from the effective field H e ff(r,f) = 
— (1 /M S )6F [m] /<5m(r, t), which is given by the magnetic 
free energy F [m]. The second term, which is param¬ 
eterized by the second-rank Gilbert damping tensor a, 
describes the magnetization dissipation, whereas the last 
term represents the current-induced torque. 

In systems with negligible intrinsic spin-orbit coupling 
(SOC)® the current-induced torque is to lowest order in 
the magnetization gradients and to linear order in the 
applied current density J given by 3 

T STT = - (i _/3m X ) (v s • V) m , (2) 

where v s is proportional to U and the spin polarization 
of the current. Importantly, we note that no effect of 
the crystal symmetry is reflected in Eq. ([ 2 ]). The torque 
is fully rotationally symmetric under separate rotations 
of the magnetization (the spin space) and spatial coordi¬ 
nates (the real space). Furthermore, the torque vanishes 
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FIG. 1: (Color online), (a)-(b) In metallic ferromagnets with 
negligible SOC, the magnetic texture is required to have a 
coupling between the itinerant quasi-particles and the mag¬ 
netization. (a) The current-driven domain wall motion is re¬ 
ciprocal to (b) the spin-motive force induced by a moving 
domain wall, (c)-(d) In ferromagnets with broken spatial 
inversion symmetry, SOC mediates a coupling between the 
magnetization and the itinerant quasi-particles, even without 
a magnetic textures, (c) The SOT on a uniform magneti¬ 
zation is reciprocal to (d) charge pumping by homogeneous 
magnetization precession. 


in a homogeneous ferromagnet because the gradients of 
the magnetization are zero in this case. The first part 
of Eq. ([ 2 ]) represents a reactive STT because it preserves 
the time reversal symmetry of the LLGS equation 0 
The second term, which is parameterized by breaks 
the time reversal symmetry of Eq. ([I]) and therefore de¬ 
scribes a dissipative STT. Its magnitude is governed by 
the spin-flip scattering rate induced by extrinsic SOC and 
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magnetic impurities. 

Recent experimental 518 and theoretical works 19 27 
have shown that the torque in Eq. © does not de¬ 
scribe the current-driven dynamics of metallic ferromag- 
nets with large intrinsic SOC. In asymmetric systems, the 
SOC combined with an external electric field produces an 
out-of-equilibrium spin density, which yields a magneti¬ 
zation torque via the exchange interaction. This torque, 
which originates solely from the SOC, is referred to as 
a spin-orbit torques (SOTs). Unlike the conventional 
STT in Eq. ©> SOTs can induce magnetization dynam¬ 
ics even in the absence of magnetization gradients, and 
their particular forms are determined by the crystal sym¬ 
metry. In Ref. l28i a general phenomenology of current- 
induced torques in systems with arbitrarily strong SOC 
was developed. Based on Neumanns principle,^ which is 
a fundamental symmetry principle of crystal physics, the 
phenomenological theory outlines a systematic scheme to 
derive all possible torque terms that are allowed by sym¬ 
metry. The formalism shows that the SOC introduces 
several new torques in addition to the reactive and dis¬ 
sipative STTs in Eq. ([2]). Additionally, the reactive and 
dissipative STTs become renormalized by the SOC, and 
the simplest extension of Eq. (|2j) to account for the effects 
of the SOC isP 

r = - 7 m x H c , (3) 

where the current-induced field H c = r/J' is determined 
by the second-rank tensor rj = [rjij], which is given by 

Vij = A ijl m k + Pijkidkmi + Pijkin.m k dim n . (4) 

Here, the tensorial forms of A^j., 0 ijkh and P ljk in 
are governed by the crystal symmetry of the system. The 
tensors A-^ and A-^. describe the reactive and dissipative 
SOTs, respectively, which are present only in systems 
with broken spatial inversion symmetry, whereas the two 
tensors Pijkin and /3 ijki represent generalizations of the 
reactive and dissipative STTs in Eq. ([ 2 ]). Here and below, 
a summation of repeated indices is implied. 

Reciprocal processes are well known in several fields 
of physics. The most familiar example may be the ther¬ 
moelectric effect, where a thermal gradient induces an 
electric current, or vice versa, an electric field produces 
a flow of heat. The linear response coefficients of recip¬ 
rocal processes are connected via the Onsager reciprocal 
relations.® Thus, by measuring one of these two pro¬ 
cesses, one also gains information about the strength and 
symmetry of the reciprocal phenomenon. 

In the reciprocal effect of the current-induced torque, 
charge currents are induced by time variations in the 
magnetization, and the forces that induce such currents 
are referred to as spin-motive forces (Fig. [l]). In the ab¬ 
sence of a large intrinsic SOC and to the lowest order in 
the magnetization precession frequency and magnetiza¬ 
tion gradients, the pumped charge current is 3136 

Ji = Pm ■ ( d t m x d,m) + /3<9 f m ■ d,jn. (5) 


Eq. 0 is reciprocal to the STT in Eq. © .(Fig- 0 -b). 
The first term, which is parameterized by P , is the re¬ 
ciprocal process of the reactive STT, whereas the sec¬ 
ond term ~ (3 is reciprocal to the dissipative STT. The 
charge-pumping effect of Eq. 0 has been experimentally 
detected in a ferromagnetic nanostrip, where an electric 
signal induced by moving domain walls was observed.® 

The general form of the current-induced torque in 
Eqs. 0-0 indicates that strong intrinsic SOC will also 
increased the complexity in charge pumping. For exam¬ 
ple, the SOTs proportional to A^ and A^ imply via the 
reciprocity relations that a homogeneous precessing mag¬ 
netization pumps a charge current in asymmetric ferro- 
magnets (Fig. 0-d). This effect, which arises solely from 
the SOC, was recently observed in strained (Ga,Mn)As 3 ' 
and has been studied in several theoretical works.®®® 
However, the effects of the SOC on the charge pumping 
in ferromagnets with an arbitrary magnetic texture (e.g., 
magnetic domain walls, skyrmions or vortices) and the 
relationship between different charge pumping contribu¬ 
tions and the current-induced torques remain unknown. 

In this paper, we extend the theory of Ref. [3S] and de¬ 
velop a general phenomenology of the coupled dynamics 
between itinerant quasi-particles and magnetization in 
metallic ferromagnets. The formalism is valid for intrin¬ 
sic SOC of any strength. Our phenomenology is based on 
Neumann’s principle, which is used to derive the possi¬ 
ble forms of charge pumping and current-induced torques 
allowed by the symmetry of the system. The Onsager 
reciprocal relations are then used to link the different 
charge pumping and torque terms. We apply the for¬ 
malism to two classes of systems: i) isotropic ferromag¬ 
nets with non-uniform magnetization and ii) homoge¬ 
neous ferromagnets with Rashba and Dresselhaus SOC. 
We show that the conventional LLGS phenomenology 
(i.e., Eqs. @ and 0 ) is obtained as the non-relativistic 
limit for isotropic systems. In these systems, the SOC in¬ 
troduces two new dissipative torques/pumping terms and 
four new reactive torques/pumping terms in addition to 
the conventional terms: to the best of our knowledge, sev¬ 
eral of these additional terms have not been previously 
reported. Importantly, these new terms lead to different 
dynamics for Bloch and Neel domain walls. For Rashba 
and Dresselhaus SOC, we derive the reactive and dissi¬ 
pative SOTs and charge pumping, and we show that the 
results are consistent with recent experimental observa¬ 
tions in strained (Ga,Mn)As epilayers. 

This paper is organized as follows. First, we introduce 
Neumann’s principle and the Onsager reciprocal relations 
in Sec. [0 before deriving the general phenomenology of 
current-induced torques and charge pumping in Sec. |III| 
The formalism is applied to isotropic and asymmetric fer¬ 
romagnets in Sec. m and Sec. [Vj A summary is provided 
in Sec. ED 
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II. FUNDAMENTAL SYMMETRY PRINCIPLES 

Our phenomenology is based on two fundamental sym¬ 
metry principles in condensed matter physics: Neu¬ 
mann’s principle and Onsager’s reciprocity relations. 
The present section gives a brief introduction to the gen¬ 
eral theory of these principles. We refer the reader to 
Refs. and [5U1 for a more comprehensive treatment. 


In the present work, we consider the linear response 
regime and treat the dynamics in the local approxima¬ 
tion, where all linear response coefficients are assumed to 
be local in space and time, i.e., the flux J(r,f) depends 
only on the value of the force f(r, t) at the space-time 
position (r,f). Then, the fluxes are related to the forces 
by the equations 

J i (Y,t)=L ij (Y)f j (Y,t ), ( 9 ) 


A. Transformation properties of tensors 

In the following discussion, it is convenient to introduce 
some basic concepts regarding tensors and vectors. 

Let {1Z £ Q} denote a set of 3 x 3 matrices, which 
represent the crystallographic point group Q. The sym¬ 
metry operation 1Z can be either be a proper rotation 
with a determinant \1Z\ = 1 or an improper rotation, 
which changes the orientation of the coordinate system 
and therefore has a determinant [R\ = — 1 . 

Vectors can be classified into two groups according to 
how they transform under a coordinate transformation: 
r = IZy. A polar vector u transforms as u = 1Zvl , 
whereas a pseudo vector transforms as u = \TZ\lZu. Im¬ 
portantly, polar vectors change sign under the spatial in¬ 
version operation TZ t j = —6ij ( Sij : the Kronecker delta), 
whereas pseudo vectors remain unchanged. Examples 
of pseudo and polar vectors are the magnetization and 
charge current, respectively, both of which are central to 
our discussion. 

Vectors are tensors of rank one. The classification of 
vectors according to their transformation properties can 
be generalized to tensors of any rank. A polar tensor 
Tij...k transforms as 

Tij...k = IZialZjp ■ ■ ■ 7(6) 

under a coordinate transformation r = 7Zy, whereas an 
axial tensor (or pseudo tensor) transforms as 

Tij...k = \R\1ZicJZjp ■ ■ ■ IZkjTai 3... 7 . (7) 

Next, we introduce the concepts of thermodynamic 
forces and fluxes, which are required in the formulation 
of Onsager’s reciprocity relations. 


B. Thermodynamic forces and fluxes 

Consider a system described by a set of thermody¬ 
namic variables {ft}. Examples of such variables include 
the charge density, spin density, and magnetization. Let 
fi denote the thermodynamic force that induces a flux 
Ji in the quantity ft. In general, the fluxes and forces 
are not uniquely defined, but for a system with uniform 
temperature, they should be properly normalized to yield 
the entropy generation S via the relation 

TS = JJi. ( 8 ) 


where the linear response coefficients Lij depend on the 
equilibrium value of the variables {ft}. The local ap¬ 
proximation implies that the magnetization slowly varies 
over typical transport length scales such as the mean free 
path, spin-flip diffusion length, and transverse-spin deco¬ 
herence length. 


C. Neumann’s principle 


The crystal symmetry of the system reduces the num¬ 
ber of independent coefficients . Formally, this reduc¬ 
tion is expressed by Neumann’s principle, which states 
that ’’any type of symmetry that is exhibited, by the point 
group of the crystal is possessed by every physical propeHy 
of the crystal”. Neumann’s principle implies that the 
response coefficients satisfy the symmetry relations 

Lij( h,m) = \n\^n ia n j0 L aP (u,M), ( 10 ) 

where Q = 0 = 1 ) if ft is a component of a vector 

(a pseudo vector) and H represents an external magnetic 
field. 


D. Onsager reciprocal relations 

The off-diagonal terms of the response matrix L = 
[Lif describe how a thermodynamic force f 3 of a vari¬ 
able qj induces a response in another thermodynamic 
variable ft. The reciprocal phenomenon of this effect is 
the process by which the thermodynamic force fi induces 
a response in the quantity qj. The linear response coeffi¬ 
cients of the reciprocal phenomena are connected via the 
Onsager reciprocal relations, which imply that 


L y (H,M) = eiejLjii- H, -M). (11) 


Here, e* = 1 (ft = —1) if qi is even (odd) under time 
reversal. 


The symmetry relationships implied by Eqs. (10 1 and 


(111 considerably simplify the possible allowed forms of 
the response matrix L and the number of independent 
matrix elements. 
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TABLE I: Fluxes and forces 

in metallic ferromagnets. 


Magnetic System 

Charge System 

Flux 

rh 

J 

Force 

fm = M s Hoff 

E 

Ci 

1 

0 

ti 

-1 

1 


III. SYMMETRY RELATIONS OF CHARGE 
PUMPING AND CURRENT-INDUCED 
TORQUES 


Eq. (13): 


L mm (m, Vm) 


— [l-0a] -1 0. 
M s 1 J 


The 3x3 matrix L pp is the conductivity tensor 

L pp (m. Vm) = (T (m, Vm), 


(16) 


(17) 


whereas the 3x3 matrices L pm and L mp describe 
the charge pumping and current-induced torque, respec¬ 
tively. 

Based on the general theory presented in Sec. [TIJ we 
derive the symmetry properties of these four submatrices. 
Under a symmetry transformation 7 Z G Q, Neumann’s 
principle implies that 


We now use Neumann’s principle to derive the possi¬ 
ble expressions for the charge pumping and the current- 
induced torque. Then, the Onsager reciprocal relations 
are used to derive relationships between these expres¬ 
sions. 

To use the Onsager reciprocal relations, the LLGS 
equation ([I]) should be rewritten such that the time 
derivative of the magnetization only appears on the left- 
hand side of the equation. We refer to this as the Landau- 
Lifshitz (LL) form of the equation. To facilitate this re¬ 
formulation, it is convenient to introduce the matrix 

Oij — tikjTYlki ( 12 ) 

where end is the Levi-Civita tensor. Then, the LLG equa¬ 
tion ([Tj) in the LL form becomes 

m = -7 [1 - Oa] _1 OH c g + r LL , (13) 
where the torque in the LL formulation is 

t ll = [1 — Oq]” 1 t. (14) 

First, we determine the different fluxes and forces. 
We consider a metallic ferromagnet at constant tempera¬ 
ture, where the magnetization is coupled to the itinerant 
charge carriers. The magnetic system is described by the 
magnetization m, whereas the charge system is charac¬ 
terized by the charge density p. The associated fluxes 
of these quantities are the time derivative of the magne¬ 
tization rh and the charge current density fj, and the 
corresponding forces that induce the fluxes are the effec¬ 
tive field M s H c ff and the electric field E. These fluxes 
and forces are summarized in Tabl e [!} where the values of 
the parameters Q and e* in Eqs. ( flO| ) are also pro¬ 
vided. In the linear response regime, these fluxes and 
forces are related via 


L mm (m, Vrrij = 7?.L mm (m, Vm) U T , (18) 

L pp (rh,Vrh) = KL PP {m,Vm)K T , (19) 

(rh,Vm) = |7^|7?.L pm (m, Vm) 1Z T , (20) 
(rh,Vrh) = \K\KL mp (m, Vm) 7L T , (21) 


L pm 


j -jmp 


whereas the Onsager reciprocal relations yield the follow¬ 
ing relationships for the off-diagonal matrix elements 


L™ m ( m,Vm) = L% m (- m,-Vm), (22) 

iij( m > Vm ) = Lji (— m , — Vm) . (23) 

L“ p ( m ,V m ) = -L”T(-m,-Vm). (24) 

The current-induced torque and charge pumping are 
usually expressed in terms of the applied charge current 
density and the time derivative of the magnetization. To 
linear order in the magnetization precession frequency 
and the out-of-equilibrium charge current density, the 
most general forms for the charge pumping and torque 


in the local approximation are 


J = fym, 

(25) 

r = - 7 m x H c , 


H c = 77 J. 

(26) 


Here, fj = [fjij] and 17 = [17^} are second-rank tensors that 
act to the right on rh and , respectively. The tensors 
depend on the equilibrium value of the magnetization and 
its gradients, i.e., fj = fj (m, Vm) and 17 = 17 (m, Vm). 
H c is the effective field produced by the current density 
and contains both dissipative and reactive contributions 
to the torque. Note that Eq. (26) refers to the current- 
induced torque in the LLGS equation ([!]). In the LL 
formulation the torque becomes 



The 3x3 matrix L mm describes the magnetization dy¬ 
namics driven by the effective held and is determined by 


r LL = -7 [1 - Oo :] -1 OH c = M s L mm r]J. (27) 

To derive the symmetry properties of the two tensors fj 
and r/, we consider an out-of-equilibrium current density 
caused by an external electric held E and magnetization 
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dynamics driven by an effective magnetic field H e g. Us¬ 
ing J = L PP E and rh = L mm f m in Eqs. ((251) and ((271), 


tensors fj and rj must transform as 


we find that fj and rj are related to the response matrices 
L pm and L mp by the equations 


L pm = fjh mm , 

L mp = M s U nm rjL pp . 


(28) 

(29) 


»7 ^m, VmJ = \lZ\TZfj (m, Vm) 1Z T , (30) 

rj (m, Vm) = \TZ\1Zrj (m, Vm) 1Z 1 . (31) 


In Eqs. the matrices L mm and L pp transform Ec l s - determine the tensorial forms of fj and rj. 


according to Eqs. (18)-(19) under a symmetry transfor¬ 


mation 1Z. Thus, for L pm and L mp to satisfy the sym¬ 
metry relationships implied by Eqs. (20) and (21), the 


In addition to the symmetry requirements imposed by 
Neumann’s principle, the Onsager reciprocal relations in 
Eq. (24) yield the equation 


M s L mm (m, Vm) rj (m, Vm) L pp (m, Vm) = -lT m (-m,-Vm) T fj (-m, - Vm) T . 


(32) 


Using Eq. (22), we identify the following reciprocal rela¬ 
tionship between fj and rj: 


M s tj (m, Vm) er (m, Vm) = — fj (— m, — Vm) . (33) 


Here, we use the conventional expression er for the con¬ 
ductivity tensor instead of L pp . 

Because the tensor M s ijcr is directly connected to the 
charge pumping tensor fj via the reciprocity relations, it 
is useful to introduce the notation 


fj (m, Vm) = M s ij (m, Vm) er (m, Vm). 


(34) 


Eq. (p6| implies that fj describes the effective field in- 

(35) 


duced by the external electric field 


1 


H c = - fj E. 

c M„ 1 


The direction of the applied electric field is typically 
known in an experiments, whereas the local current den¬ 
sity is determined by local magnetoresistance effects. 
Therefore, it is often preferable in systems with large 
SOC to express the current-induced torque in terms of 
the electric field (using the tensor fj). 

Using the transformation properties of the conductiv¬ 
ity tensor and ij , one finds that fj satisfies the transfor¬ 
mation rule 


fj ^m, Vmj = \TZ\lZfj (m, Vm) 1Z T , 


(36) 


Next, we will use the relations (301, (31), and (36) with 
Eqs. (|25|)-(26) to develop a systematic scheme to derive 


all possible torque and charge pumping terms. 


A. Phenomenological expansion 


In general, several tensors satisfy the symmetry re¬ 
quirements imposed by Eqs. (30), (31) and (36). This 


situation is analogous to the case of finding the mag¬ 
netic free energy that describes the magnetocrystalline 
anisotropy energy. In that case, one has a free en¬ 
ergy functional that satisfies a set of symmetry rela¬ 
tions F [m] = F[\TZ\TZm] generated by the point group 
{1Z £ Q}. It is common to approximate the free 
energy by performing an analytic expansion F [m] ss 

K^rriimj + K^ kl mimjm k mi +... in powers of m, where 
each term in the expansion satisfies the correct symmetry. 
Usually, the first harmonics in this expansion capture the 
magnetocrystalline anisotropy observed in experiments. 
The series expansion rapidly converges because higher- 
order terms have a higher degree of anisotropy and thus 
tend to average out due to thermal fluctuations. Addi¬ 
tionally, the higher-order terms are believed to be of a 
higher order in the SOC. 

With this motivation, we perform a similar analytic 
expansion of the tensors fj , fj and rj: 

fjij = + A [f k m k + j3ijkid k mi + Pij k i n m k dim n , 

fjij = A^ } + A [f k m k + hjkid k m + Pijkinm k dim n , 

rjij = + A \f k m k + Pijkid k mi + Pij k i n m k dim n . 

(37) 


In Eq. (37), we keep only the first harmonics, which rep¬ 
resent generalizations of the torque and charge pumping 
expressions in Eqs. ([2]) and ([5]). The expansion can be 
continued to arbitrary high orders in m to fit the experi¬ 
mental observations for a particular system. Such higher- 
order harmonics represent the current-induced torques 
and charge pumping with higher degrees of anisotropy. 


The Onsager reciprocal relations (33) imply the follow- 
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ing relationships among the tensors in Eq. (37) 


A (r) 

*7 

= —A (r) 

V ’ 

(38) 

A (d) 

9-5 

'<1 

II 

(39) 

fiijkl 

— fijikh 

(40) 

Pijkin 

Pjikln • 

(41) 


Eqs. (38)-(41) connect the reciprocal terms in the phe¬ 


nomenological expansions of the charge pumping and 
current-induced torque. For example, A^ represents the 
charge pumping process that is reciprocal to the torque 
given by the tensor A . 

In general, the conductivity tensor also depends on the 
local direction of the magnetization and its gradients. 
To capture this dependency, one can perform a similar 
phenomenological expansion of the conductivity in pow¬ 
ers of m. An expansion of the diagonal tensor elements 
describes the anisotropic magnetoresistance (AMR) ef¬ 
fect, whereas an expansion of the off-diagonal elements 
describes the anomalous Hall effect (AHE). By relating 
terms of equal order in m in the two expressions M s rjcr 
and fj, one can derive relationships among different ten¬ 
sor elements in the expansions of r) and fj. 

The forms of the tensors in Eq. (37) are determined 
by Eqs. (30)-(31) and (36). Because the tensors fj , fj 


and r) transform similarly under a coordinate transfor¬ 
mation, the terms in their phenomenological expansions 
also satisfy identical symmetry relations. Therefore, it is 
sufficient to write down the relationships satisfied by the 
tensors in the expansion for the current-induced torque: 


A (r) = 

a 12 = 


a%, 


ijk 

fiijkl 

Pijkln 


33 kk ijk 

'k'ii'kjj , 'kkk l 'R'U’fii'j'k'l' i 

'kj j' 'R'kk' 'k-ll' 'k' 


,P, , 

nn i j k L n 


(42) 

(43) 

(44) 
• (45) 


a (r) 

From these transformation rules, we observe that A\- 
and Pijkin are axial tensors, whereas Aand fiijki are 

polar tensors. Importantly, we note that A^ = A^j. = 0 
in systems that are invariant under spatial inversion 
(i.e., under lZ,j = —Sij). Thus, as it is established, 
charge pumping and current-induced torques are ab¬ 
sent in homogeneous ferromagnets with inversion symme- 
tryPl22l223i 2 ' Only magnetization gradients can mediate 
a coupling between the itinerant charge carriers and the 
magnetization in such systems. The tensors AA 

Aq\ MjL Aand A.\f k play a role exclusively for crys¬ 
tals described by a non-centrosymmetric point group. 

Next, we apply the general theory developed in this 
section to two classes of systems: isotropic systems and 
asymmetric systems with Rashba and Dresselhaus SOC. 
An example of an isotropic system with SOC is the disor¬ 
dered alloy Permalloy, whereas the ferromagnetic semi¬ 
conductor (Ga,Mn)As is an example of a system with 
both large Rashba and Dresselhaus SOC. 


To circumvent the complications introduced by mag¬ 
netoresistance effects (as mentioned above), we will ex¬ 
press the torques in terms of the electric field instead of 
the charge current density. 


IV. ISOTROPIC SYSTEMS 


Isotropic systems are invariant under any proper or 
improper rotation. Thus, spatial inversion is a symmetry 
operation, which implies that 


A (r ) - A (d) - A (r) - A (d) - A (r) - A (d) - 0 

ly ij ~ ly ijk ly ijk ~~ ly ij ijk u ' 


(46) 


Therefore, the coupled dynamics of the itinerant quasi¬ 
particles and the magnetization is governed by the last 
two terms of the expansions in Eq. (37) (parameterized 
by the /3 and P tensors). 


A. Invariant tensors of isotropic systems 


It is well known that any isotropic tensor of rank two 
is proportional to the Kronecker delta 5ij , whereas all 
isotropic third-rank tensors are proportional the Levi- 
Civita tensor e l3 k- Furthermore, it has been shown that 
all isotropic tensors of even rank are linear combinations 
of products of Kronecker deltas, whereas isotropic ten¬ 
sors of odd rank are linear combinations of terms formed 
from products of the Kronecker delta and the Levi-Civita 

tensor.^ 


Because the (3 and P tensors in Eq. (37) are of rank 
four and five, respectively, we concentrate on these two 
types of tensors. In general, an isotropic tensor of 
rank four has three independent tensor coefficients and 
can be written a^^ 


1§h = T^SijSki + T^StkSji + T^SuSjk- (47) 


An isotropic tensor T^ kln of rank five has six independent 
coefficients ^ 


1ijkln ~ ^ e ijk<3ln + e ijl$kn + Cijn^kl + 

2"4 ^CiklSjn + 2~5 ^ikn^lj + Tq ^iln^jk- ( 48 ) 


Consequently, fiijki and /3ijki have tensorial forms given 
by Eq. (47), whereas the tensorial forms of Pijkin and 
Pijkin are determined by Eq. (48). 

We will later on apply the following useful identity ^ 


tijk^mp tjkm^ip T ^-kmi^jp ^mijOkp 


j&kp — 0 , 


(49) 


which can be used to decompose an isotropic tensor of 
rank five into a linear combination of other terms. 
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B. Expressions for torque and charge pumping 

In the following discussion, it is convenient to sepa¬ 
rate the charge pumping and torque terms into dissipa¬ 
tive and reactive contributions. The reactive terms arise 
from Pijkin and Pijkim whereas the dissipative terms are 
determined by the tensors $%jki and ftijki- 

Using Eq. (47), we find the dissipative current-induced 


held and charge pumping jf*' ' 

— /3im ■ dim + ■ V?rq + /3 3 to,;V ■ m, (50) 

1 


= 


L 


(i-\ E ■ Vmj + /3 2 E • d t m + ft 3 E t V ■ m 


The reactive held He ■* and charge pumping ^ each 
have six independent terms. However, because of the 
normalization condition m • m = 1, one of these terms 
vanishes and the resulting expressions for the reactive 
held and charge pumping become 


j, m = 


h { :) = 


Pim • (m x 9im) + P 2 (m x m) ■ (V • m) 

+ P 3 (m x (m • V) m) i + f 4 m • (m x V) m,; 

— P 5 rnjrii • (V x m), (51) 

^-[Pi (m x (E • V) m) i + P 2 (V • m) (E x m). 

+ P 3 (E x (m ■ V) m)- + P 4 E ■ (m x V)^ m 
+ P 5 (E • m) (V x m).]. (52) 


The Onsager reciprocal relations (40) - (41) imply that 


the phenomenological parameters in Eq. (50) - (52) are 
related via the equations 


on the spin space and the real space , i.e., acting on both 
the spatial vectors (such as VU and p) and the spin vec¬ 
tor <x. However, without SOC, the spin space and real 
space are no longer coupled and can therefore be sepa¬ 
rately rotated. This decoupling increases the symmetry 
of the system: an isotropic system without SOC is invari¬ 
ant under separate rotations of the spin space and the real 
space. 

The expressions in Sec. 


IV B are invariant under si¬ 


multaneous rotations of spatial vectors and spin vectors, 
i.e., they apply to isotropic systems with SOC. If we also 
require that the expressions are invariant under separate 
rotations of the spin space and real space, we hnd that 
only the dissipative contributions proportional to fti , j3\ 
and the reactive contributions proportional to P \, Pi sat¬ 
isfy the correct symmetry. Without SOC, the conductiv¬ 
ity tensor becomes cy,- = aSij and 17 = M s ar]. Thus, for 
an isotropic system without SOC, the charge pumping 
and torque are given by 


Ji = Pim • (m x <9,m) + j3\va • c^m, 


(56) 




Eq. (56) is identical to the charge pumping in Eq. ([ 5 ]) , 
whereas Eq. ([57|) reduces to the STT in Eq. 0 when 


7 1 P 1 J/M s a = — v s and j/3iJ/M s a = -/3v s . 

From the above analysis, we conclude that the conven¬ 
tional LLGS equation 0 and charge pumping 0 are ob¬ 
tained as a particular limiting case for our phenomenol¬ 
ogy: the non-relativistic limit of an isotropic system. 


D. The relativistic corrections 


h=k i e {1,2,3}, (53) 

Pi=Pi, ie{l,...,5}. (54) 

Now, we will demonstrate that the conventional LLGS 
equation ([l]) and charge pumping ([ 5 ]) represent a partic¬ 
ular limiting case for isotropic systems. 


C. Non-relativistic limit: conventional LLGS 

To illustrate how the system is affected by the SOC, 
we consider the SOC Hamiltonian 


iPn = 


4to 2 c 2 


(W x p) ■ a. 


(55) 


Here, V(r) is the crystal potential, p is the momentum 
operator, <7 is the vector of Pauli matrices, m is the elec¬ 
tron mass, and c is the speed of light. Importantly, H so 
vanishes in the non-relativistic limit c —> 00 . The Hamil¬ 
tonian ( |55| links the spin space to the symmetry of the 
underlying crystal structure. For an isotropic system, 
Eq. (55) is (after disorder averaging) invariant under any 


proper or improper rotation TZ that acts simultaneously 


Our phenomenology shows that the SOC introduces 
two new dissipative terms and four new reactive terms in 
isotropic systems such as Permalloy. Thus far, none of 
these terms have been inferred from experiments. How¬ 
ever, if these terms are comparable in magnitude to the 
conventional, non-relativistic terms, they have important 
consequences for the current-driven dynamics of domain 
walls. Because V • m = 0 (V • m 7 ^ 0) for Bloch (Neel) 
walls, the terms proportional to P 2 and ft 3 only con¬ 
tribute to the dynamics of Neel walls, which is also the 
case for the other new terms; the terms contribute differ¬ 
ently for different types of domain walls. Therefore, even 
in isotropic systems, Bloch and Neel walls are expected 
to exhibit different current-driven drift velocities. 

A way to probe some of the relativistic correc¬ 
tions is to investigate the current-driven domain wall 
drift velocity for electric fields applied parallel to the 
wall. In this case, only the relativistic corrections to 
the reactive STT is able to drive domain wall mo¬ 
tion. Let us consider Bloch and Neel walls of the 
forms m (z,t) = [sisech(5), s 2 tanh(z), 0 ] and m (z,t) = 
[0, sisech(z), s 2 tanh(z)], respectively. Here, z = (z— 
rdw )/A where r^w is the domain wall position along the 


















z=[001] 


b 


z axis, A is the domain wall width, and Si G {—1,1}, 
S 2 £ {—1,1}. Using the collective coordinate descrip¬ 
tion outlined in Ref. |43J we find for an electric field 
E = [E x . E y , 0] the following stationary drift velocities 
for Bloch and Neel walls 


Bloch : 

- 7rP 4 „ 

t dw — '■’i/i Py> 

4 a 

(58) 


7T P2 

(59) 

Neel : 

^ dw — ^1 Px • 

4a 


In a real experiment, it is difficult to produce an electric 
field that is perfectly parallel to the wall. Therefore, in 
practice, there will most likely also be a small E z com¬ 
ponent. This E z component can induce a drift velocity 
via the conventional dissipative STT, i.e., f,i w ~ /3/a. 
However, Eqs. (58)-(591 implies that a reversed domain 
wall, m — > — m, moves in the opposite direction. This 
differs from the drift velocity induced by the dissipative 
STT, which is unchanged under m — ► — m. Thus, if mea¬ 
surements of up-down and down-up domain walls yield 
opposite drift velocities, this is a clear indicator that the 
relativistic reactive STT is the driving mechanism behind 
the current-induced velocity. 

Some of the relativistic corrections in Eqs. (50)-(52) 
can be attributed to the spin Hall effect (SHE). Ref' 1441 
showed that the SHE produces reactive and dissipative 
torques of the form 


T she = x m) • Vm, (60) 

T^he = ct H m X [(E x m) • Vm], (61) 


Here, an is the Hall angle. Using the identity in Eq. (491, 
the dissipative SHE torque can be decomposed into 


T she,i = -P<xh[ E • (rn x V^m - (m x E)*(V • m)]. 


The first term is of the same form as the /3 2 torque in 
Eq. (50), while the second term has the form of the /3 3 
torque. Taking the cross product with m, one finds from 
the above expression that reactive torque can be 


written as a linear combination of the P 2 and P 4 terms in 
Eq. ( [52| ). The reciprocal processes of these SHE torques 
were calculated microscopically in Ref. 1451 

We believe an important task for future experiments 
is to investigate the magnitude of the relativistic torques 
in isotropic ferromagnets such as Permalloy. Addition¬ 
ally, improved microscopic theories are required to gain 
a deeper understanding of the relativistic corrections to 
the current-driven magnetization dynamics and to deter¬ 
mine the magnitude of the new terms. 


a 



[ 110 ] 
X X 


{[ 110 ] 




FIG. 2: (Color online), (a) A (Ga,Mn)As epilayer grown on a 
GaAs substrate along the crystallographic direction [001]. (b) 
Stereographic projection of the point group C 2 V , which illus¬ 
trates the symmetry of strained (Ga,Mn)As. (c) Vector plot 
of the reactive SOT fields induced by Rashba (red arrows) 
and Dresselhaus (blue arrows) SOC for different directions 
of the applied electric field, (d) Vector plot of the in-plane 
current pumping due to magnetization precession for differ¬ 
ent directions of the equilibrium magnetization. The charge 
pumping via Rashba (Dresselhaus) SOC is illustrated by red 
(blue) arrows. The reactive and dissipative contributions to 
the pumped in-plane current have identical symmetries. 


grown on GaAs along the crystallographic axis [001] 
(Fig. [§1 -b). In this case, the two-fold symmetry axis is 
perpendicular to the epilayer (i.e., along [ 001 ]), and the 
two reflection planes are defined by the axes [ 110 ] and 
[110]. Herein, we define the reference frame such that 
the two-fold axis is along z, whereas the x and y axes 
define the two reflection planes. For (Ga,Mn)As, this 
reference frame is related to the crystallographic axes 
via x = [110], y = [110], and z = [001]. We consider 
ferromagnets with a uniform magnetization and model 
the coupled dynamics of the magnetization and itiner¬ 
ant charge carriers using the four second- and third-rank 
tensors A^, A^, A< d ), and A {d) in Eq. ((37]). Note that 
these terms originate solely from the SOC and vanish 
in systems with spatial inversion symmetry. Thus, A^ 
and A^) parameterize an SOT, whereas Ay> and A ^ 
describe charge pumping caused by a direct conversion 
of zero wavevector magnons into an electric current via 
the SOC, which we previously termed magnonic charge 
pumping . 37 


A. Invariant tensors of the C^v point group 


V. RASHBA AND DRESSELHAUS SOC 

Next, we consider asymmetric ferromagnets with 
Rashba and Dresselhaus SOC. These systems are de¬ 
scribed by the crystallographic point group C^v- An im¬ 
portant example of such a system is strained (Ga,Mn)As 


( 2 ) 

The invariant axial second-rank tensors of the 

l J 

point group C 21 , are described by two independent tensor 
coefficients. For the reference frame defined above, the 
two non-vanishing tensor coefficients ard^H 

/ri(2) rp( 2) 

xy 5 yx ’ 


( 62 ) 
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(3) 

Polar third-rank tensors that are invariant under C 2 „ 
are described by the following seven independent tensor 
coefficient^^! 


Ji(3) rp( 3) rp(3) rp(3) rp( 3) rp( 3) rp( 3) 

zzz"> M xxzi ^ zxx") ^xzx") - L yyz'> ^ zyyi ^yzy’ \ yju / 


Thus, charge pumping and SOT are described by nine 
independent tensor coefficients. The tensors AM and 
are characterized by the two elements in Eq. ( [62] ) , 
whereas AM and AM are determined by the seven coef¬ 
ficients of Eq. (63). 


B. Expressions for torque and charge pumping 

We separate the various terms into reactive and dissi¬ 
pative contributions. The reactive charge pumping and 
SOT field are determined by the second-rank tensors AM 
and AM, respectively, which lead to the expressions 


j} r) = 6 ix A^m y + S iy A^rr 


M s H C!i 


ly ^ yx m x , (64) 
M _ V Ur) t? I r. UU TP 

— U IX 1Y T.1I V I u IV 1y ht. J-'X- 


*-xy 


Jr t'V ± '-yx 


(65) 


The dissipative charge pumping and SOT are given by 
the third-rank tensors AM and AM, respectively. From 
Eq. (63), we find the pumped current and the effective 
held: 


-7.M _ 


AM 2 m x m 2 + AM,to 2 to x 


viy 


AM 


M S H { C J = 4 


y yy Z m y m z + AM \rn z m y 


u iy 


^xxz^x m z + ^xzx^z m x 

^yyz-Ey m z S' ^yzy^z m y 


( 66 ) 


AM rh z m, + AM rh x m x + AM rf 


m v m , 


(67) 


+ 


A^J z E z m z + AM ,E x m x + AM,£, 


zxx J -‘x" u x 


s -zyy J ^y m y 


The reciprocal relationships between the phenomeno¬ 
logical pumping and torque parameters are determined 
by Eqs. (38) - (39). 


As an application of Eqs. ]64fo - (67), we focus on the 


specific system of (Ga,Mn)As, for which the SOT and 
charge pumping have been experimentally observed. 3 ' 


C. Example: (Ga,Mn)As epilayers 

We consider a strained (Ga,Mn)As epilayer grown on 
a GaAs [001] substrate. At equilibrium, the magnetiza¬ 


tion is assumed to align along an in-plane direction (i.e., 
perpendicular to [001]). In addition, we only consider 
in-plane electric fields and pumped charge currents: 


m = [m x m y 0], E = [E x E y 0], J = \j x J v 0 . 
We parameterize the direction of m and E by the angles 
<j> m and <j) e , respectively, such that E x = Eq cos((p e ), E y = 
Eq sin(0 e ), m x = cos(<^ m ), and m y = sin(</> m ). These 
angles are measured from the x axis, which is parallel to 
the crystallographic direction [110]. For later use, we also 
introduce the angle (f> me , which represents the relative 
angle between the magnetization and the applied electric 

field 0 rn e = 0rn 0e ■ 


Because the magnetization and applied electric field 
do not have a z component, several terms in Eqs. (64) - 
(67) vanish. The reactive charge pumping and current- 
induced field remain characterized by the same two ten¬ 
sor coefficients, as given above. However, the dissipative 
contributions become considerably simplified. The dis¬ 
sipative (in-plane) charge pumping is governed only by 
the two terms that are proportional to A[?, and A y %, 
whereas the dissipative SOT field is given by A£t® and 
Ai d yy . Because the pumping and torque parameters are 
related via the reciprocal relations, we conclude that only 
four independent parameters are required to model this 
system. It is convenient to decompose the four terms into 
Rashba-like and Dresselhaus-like contributions by intro¬ 
ducing the parameters 


Afl } = (AM - AM)/2, A£> 
A^ = (Ait + AM)/ 2 , A^ 
A^ = (AM - AM)/2, A£> 
A^ = (AM +AM)/2, AM 


(AM+AM)/2, 
(Ait - AM )/ 2 , 
(AM+AM)/2, 
(Ait-AM)/2. 


( 68 ) 


Here, AM (AM) and A^ (A^) parameterize the reac- 


VA?)— 

tive SOT (pumping) caused by Rashba and Dresselhaus 
SOC, respectively, whereas A^ (A^) and A^f (A^f) 
represent the corresponding dissipative processes. By 


describing Eqs. (64) - (67) in terms of the parameters 
in Eq. (68), we obtain the expressions 
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H r . z = 


The charge pumping of Eq. (69) was recently investi¬ 


gated and measured in Ref. [371 The reactive and dissipa¬ 
tive charge pumping with Rashba symmetry agree with 
the results derived microscopically for a Rashba model in 
Refs. 1381401 The reactive parts of the SOT in Eq. (f70j) 
have been derived microscopically in several works,- 172 
while the dissipative SOT © agrees with the results 
derived in Refs. 11 71241261 

The charge-pumping and reactive SOT-held symmetry 
of Eqs. ( [69] ) - (pflj) can be represented in a vector plot. 
Fig# -d shows the direction of the reactive and dissipa¬ 
tive pumped charge current for different directions of the 
equilibrium magnetization and the direction of the reac¬ 
tive SOT held for different directions of the applied elec¬ 
tric held. Recent experiments have verified that the SOT 
and charge pumping in (Ga,Mn)As follow the f orms of 
Eqs. ( |69| - ( [7l| ), as implied by our phenomenology.^® 17 37 

Because both reactive and dissipative processes con¬ 
tribute to the in-plane charge pumping, these processes 
can be difficult to disentangle from a measurement of 
the charge pumping. However, for the SOT held, only 
the reactive processes control the in-plane components 
of the held, whereas the out-of-plane component H c z is 
governed by the dissipative effects. 

The magnitude of H c z can be written as a function of 
two angles, <j) me and (j) e . From Eq ([71]), we hnd 


A (d) A (d) 

~^~E 0 cos(0 me ) + -jj-E 0 cos(0„ 


■20 e ). (72) 


The (f>me dependence of H CjZ is illustrated in Tab. [TT] for 
four different directions of the applied electric held and 
is consistent with the out-of-plane held observed in Ref. 

nn 

The reactive and dissipative charge pumping caused by 
Rashba (Dresselhaus) SOC are related via the Onsager 
reciprocal relations in Eqs. (38) - (391: 


A ); 1 = -A« 

A« = kf, 


A (r) - -A (r) 
1\ D ~ > 

A (d ) - A (d) 

1V D ~ 1V D * 


(73) 

(74) 


Thus, by measuring one of these two processes (e.g., the 
SOT), one also automatically gains knowledge regard¬ 
ing the strength of the reciprocal phenomenon (e.g., the 
charge pumping). This reciprocal relationship between 
magnonic charge pumping and SOT was recently demon¬ 
strated experimentally in (Ga,Mn)As.® 7 


TABLE II: Angle (/> me dependence of the Rashba and Dres¬ 
selhaus contributions to the SOT held component H CiZ . 



0e 

Rashba SOC 

Dresselhaus SOC 

E| | [100] 

7T 

4 

cos(</> me ) 

sin(</> me ) 

E|| [010] 

7r 

4 

COS (0 me ) 

- sin(0 me ) 

E||[110] 

0 

COS (0 me ) 

COS (0me) 

E| | [110] 

7T 

2 

COS (cpme) 

- COS(0me) 


VI. SUMMARY 


Several experiments have revealed an extreme com¬ 
plexity in the coupled dynamics between itinerant charge 
carriers and magnetization in systems with strong SOC. 
The theoretical formalism in this paper outlines a system¬ 
atic scheme for deriving all possible torque and charge 
pumping terms that satisfy the correct symmetry im¬ 
plied by Neumann’s principle and relates these terms 
via the Onsager reciprocal relations. We hope that the 
present formalism will play an important role in the ex¬ 
ploration of these systems and will be a useful tool for 
understanding and directing future experiments on fer- 
romagnets with strong SOC. 

We apply the formalism to two important classes 
of systems: isotropic ferromagnets with non-uniform 
magnetization and homogenous ferromagnets described 
by the crystallographic point group C^,. In isotropic 
systems, the phenomenology reduces to the conven¬ 
tional LLGS equation and spin-motive force in the non- 
relativistic limit. However, we show that the SOC in 
these systems introduces six new reactive and dissipative 
terms. Four of these new terms can be attributed to the 
SHE, whereas two of the reactive terms have not been 
reported before. Importantly, the new terms lead to dif¬ 
ferent dynamics for Bloch and Neel domain walls, i.e., the 
dynamics are dependent on the structure of the magnetic 
texture. For asymmetric ferromagnets described by the 
point group C^v, we concentrate on strained (Ga,Mn)As. 
We show that the SOTs and charge pumping derived 
from our phenomenology are consistent with recent ex¬ 
perimental observations in (Ga,Mn)As epilayers. 
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